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Abstract 

We have performed the most comprehensive predictions of the temperature fluctuations ST/T 
in the primeval isocurvature baryon models to see whether or not the models are consistent with the 
recent data on the cosmic microwave background anisotropies. More specifically, we computed the 
ST/T corresponding to the experimental set-up by the South-Pole and the Owens Valley experiments 
as well as the COBE satellite. The amplitudes of the predicted ST/T are normalized by means of 
the COBE 10° data. The resulting constraints on the models are presented on n — fib plane in the 
case of Ao = 1 — J2& (flat models) and Ao = (open models), where n is the primordial spectral 
index of entropy fluctuations and is the present baryon density parameter. Our results imply 
that the PIB models cannot be reconciled with the current observations for any reasonable set of 
cosmological parameters. 
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1. Introduction 

The discovery of the temperature fluctuations ST/T in the cosmic microwave background (CMB) 
by COBE (Smoot et al. 1992) strongly supports the gravitational instability picture for origin of 
hierarchical structures in the universe. In particular, the amplitude of ST/T places tight constraints 
on the model parameters of viable scenarios (Wright et al. 1992; Gouda & Sugiyama 1992; Sugiyama 
& Gouda 1992). In the case of the standard cold dark matter (CDM) model in which the total density 
parameter fio is unity, the dimensionless cosmological parameter Ao vanishes, and the dimensionless 
Hubble constant h (= _H"o/100km/sec/Mpc) is 0.5, for instance, the COBE ST/T on 10° scale implies 
that the amplitude of top-hat mass fluctuations averaged over r = 8/i _1 Mpc sphere should be close 
to unity. It has been pointed out, however, that such unbiased standard CDM models have serious 
troubles in explaining the galaxy-galaxy two-point correlation function, and small-scale pair-wise 
velocity dispersions of galaxies among others (Davis et al. 1985; Suto & Suginohara 1991; Ueda, 
Itoh & Suto 1993; Suto 1993). In the present paper, we will re-examine the viability of yet another 
model of structure formation, primordial isocurvature baryon (PIB) models (Peebles 1987a, b), in 
the light of the COBE (Smoot et al. 1992), the recent ACME-HEMT South-Pole (SP : Gaier et al. 
1992; Schuster et al. 1993) and Owens Valley (OVRO : Readhead et al. 1989) data. 

There are a number of previous studies considering the astrophysical consequences of, and 
the resulting constraints on the PIB models(Peebles 1987a, b; Efstathiou & Bond 1987; Efstathiou 
1988). They indicate that successful PIB models require a steep spectrum of primordial entropy 
(or baryon number) fluctuations; if the spectrum of the primordial entropy S(k) (k is the comoving 
wavenumber) is given by a simple power-law (|S"(A;)| 2 ) oc k", then n ~ —1 is preferred rather than 
n = —3 predicted by conventional inflation models (Gouda, Sasaki & Suto 1989; Suginohara & Suto 
1992; Cen, Ostriker & Peebles 1993). It turns out, however, that this contrived spectrum can be 
realized in a model of baryogenesis in the context of inflation (Yokoyama & Suto 1991). Unless a 
significant fraction of the present universe is dominated by non-baryonic dark matter which does 
not interact with radiation field, small scale excessive ST/T is inevitably produced both in adiabatic 
and isocurvature models (Gouda, Sasaki & Suto 1987, 1989). This is not problematic in PIB models 
in which reionization of the universe at high redshifts is expected to occur due to the considerable 
power of fluctuations on small-scales. This early reionization, however, would produce a secondary 
ST/T via the Doppler effect of the ionized plasma (Ostriker & Vishniac 1986; Vishniac 1987). 

Detailed examination of the CMB anisotropies incorporating the above velocity-induced tem- 
perature fluctuations in the PIB models was made by Efstathiou and Bond (1987), and Efstathiou 
(1988). In the present work, we solve the Boltzmann equation directly up to the present epoch 
taking full account of all the first-order effects without approximation. Then we separately com- 
pute the contribution from the second-order velocity terms discussed by Vishniac (1987) which is 
added in quadrature to the direct numerical results to the first order. In doing so, we properly 
evaluate the growth rate of the second-order terms in fio < 1 universes. With the above scheme, we 
compute the ST/T which should be compared with the COBE 10° scale and quadrupole data, the 
SP experiments, and the OVRO experiment. The resulting constraints on the models are presented 
on n — fib plane (in the present baryon-dominated models, fio = J2& ) in the case of Ao = 1 — J2& 
(flat models) and Ao = (open models). The former has not been studied quantitatively in the PIB 
scenarios so far. 
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2. Scenarios and parameters in PIB models 

2.1. Basic scenario 

In the present model, the primordial baryon number fluctuations are responsible for the forma- 
tion of cosmic structures at much later epochs. The spectrum shape at the initial epoch a,- when the 
fluctuations were generated is largely model-dependent and we simply adopt a single power-law form 
for the primordial entropy power spectrum (|S"(A;)| 2 ) oc k" which is a reasonable approximation to 
the proposed models (Yokoyama & Suto 1991; Sasaki & Yokoyama 1991) on scales of astrophysical 
interest. As mentioned earlier large-scale galaxy distribution is consistent with n = — 1 (Suginohara 
& Suto 1992; Cen et al. 1993) while conventional inflationary models predict n = —3. Since the 
universe is assumed to be dominated by baryons only, i.e., J2q = Silk damping (Silk 1968) would 
smear out the small-scale density fluctuations if they were adiabatic. This is why one has to resort 
the isocurvature fluctuations for a realistic reionization scenario, even if the mode seems somewhat 
contrived. The details of the reionization history are quite uncertain or at least model-dependent, 
and thus we assume that the universe never recombined for simplicity. This corresponds to the case 
of the maximum erasure of the primary ST/T . 

If the ionization rate of the hydrogen x e has been constant since the redshift z, the corresponding 
optical depth for the electron scattering is expressed as: 

2r 



2 2 — 3fio + (fioz + 3fio — 2)y/l + fioz (open models) 



n (i + z) 3 + i-n 



(2.1) 



(flat models) 



3O 

where To is the optical depth in the present horizon size (= c/Ho) and is given by 

3c (JrHn ( Yu e \ , 
t = — -x e Q b 1 — ~ 0.06x e Q b h. (2.2) 

In the above expression, c, <7y, G, and m p denote the velocity of light, the Thomson cross section, 
the gravitational constant, and the proton mass, respectively. In numerical evaluation of To, we 
adopt Yjj e = 0.23 as the abundance of helium by mass. The epoch when t(z) reaches unity is most 
important in discussing the effect of reionization. If the redshift of the epoch z c is much larger than 
l/fio, one may derive the following approximate expression for both open and flat models: 

/O \ 1/3 / 0.05 \ 2 / 3 , , 

z c ~ 30 — . (2.3) 

In what follows, all quantities defined at this epoch z c are indicated by the subscript c. For example, 
t c denotes the cosmic time corresponding to the redshift z c . Clearly almost complete ionization 
(x e = 1) is necessary for the past universe to be opaque with respect to Thomson scattering. 

Early reionization of the universe z z c would substantially smooth out the intrinsic anisotropies 
in ST/T , but in turn produce a significant amount of secondary anisotropies due to the Doppler 
and the inverse Compton effects (Zel'dovich & Sunyaev 1969) since some initial time i,- up to the 
present to: 

ST f*o / N 

~T = it 1 ' w1le(JTdt > ( 2 - 4 ) 
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ST [*o k B T e 

— = -2 / -n e a T cdt, (2.5) 

T Jti m e c l 

where ks is the Boltzmann constant, T e , m e , and n e are the temperature, mass and number density 
of electrons, and 7 is the direction of line- of- sight. 

The ST/T due to the cumulative Sunyaev-Zel'dovich effect (eq.[2.5]) from nonlinear objects at 
high redshifts was considered by Makino and Suto (1993). Although their results indicate that the 
effect potentially induces the observable amount of ST/T , the quantitative values are dependent 
on the assumption of the number density and the spatial correlation of the nonlinear objects at 
high redshifts. Therefore we focus our attention on the velocity-induced term (2.4) in what follows. 
Ostriker and Vishniac (1986) for the first time pointed out the possibility that the second-order 
velocity induced terms in eq.(2.4) might dominate the temperature fluctuations in some cosmological 
scenarios on small scale. Vishniac (1987) worked out that idea in detail and computed the ST/T in 
CDM models with early reionization. Subsequently, this method was refined and applied to the PIB 
models (Ao = 0) by Efstathiou and Bond (1987). Recently, Hu, Scott and Silk (1993) and Dodelson 
and Jubas (1993) evaluated the magnitudes of terms contributed to the ST/T up to the second- 
order, and concluded that the second-order velocity term considered by Vishniac (1987) gives the 
dominant contribution on small scales. 



2.2. Transfer function of baryons 

The evolution of isocurvature fluctuations is completely different from their adiabatic counter- 
part. For quantitative comparison of the ST/T predicted in the PIB models and the observational 
data, one has to numerically solve the linearized Boltzmann equations which describe the general 
relativistic evolution of photons, baryons, and massless neutrinos. Nevertheless it is instructive to 
summarize beautiful analytical results based on perturbation theory by Kodama and Sasaki (1986), 
who considered baryon-photon universes neglecting the contribution of neutrinos. They showed 
that the entropy fluctuations S(k) is nearly constant in isocurvature models, and that the spectrum 
of total energy density fluctuations after the matter-radiation equality epoch reduce to 

f <|S(A;)| 2 ) (k~>a eq H eq ) 

( 225 \a J nl(a H )4" [ ^' [ ^ e i e ^ 

where a is the scale factor at that epoch, and a eq and ao are those at the equality and the present 
epochs, respectively. Therefore even if the primordial entropy fluctuations (|S"(A;)| 2 ) are given by a 
scale-free form oc k n , the resulting total fluctuations at a > a eq develop a peak near a comoving 
scale for the Hubble length at a eq : 

\eak = TT ^ . r , ~ 45(O b /i 2 )- 1 Mpc. (2.7) 



a. 



,^(15/8^2) 



Since the present spectrum of baryons would correspond to eq.(2.6) at the decoupling epoch of 
baryons and photons, z c , the height of the peak is characterized by the ratio of the power spectrum 
at each side: 



2250 2 c V8V2 a #o / V 0.01 



2 



R s _* -» f " Sa^V „ 40 M! " . (2 . 8) 
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Therefore as fi^/i 2 increases, the more prominent peak forms at the smaller scale. The peak height 
becomes larger if the decoupling occurs late. 

The above features should be compared with the direct numerical computation. In Figure 1 are 
plotted the transfer functions of the baryon density fluctuations T(k) up to the present epoch a®: 

T(k) = (IMMo)l 2 ) , , 

Figure la displays the results for open models, while Figure lb for flat models. Curves in each 
panel correspond to = 0.1, 0.2, 0.4, 0.6, 0.8 and 1.0 from top to bottom, both for h = 0.5 and 
1.0. Since this spectrum plays a fundamental role in computing the velocity induced ST/T below, 
we fitted the numerical results (solid curves) by the following form (dotted curves): 

T(k) = [ak 2 exp (-bk 2 ) + exp (--jL + ] 2 , (2.10) 

with 

6000 _ 2 , 

m Mpc 2 , (2.11) 



(Sl b h 2 ) 



(^Ibh 2 ) 
220 



180 

Y^Mpc 2 (Q b h 2 > 0.15), 



2.12 



(fW) 



^Mpc 2 (Q b h 2 < 0.15), 



c = 0.01 (jVi^^Mpc -2 , (2.13) 

d = 0.02 (jVj 2 ) ' 8 Mpc -1 . (2.14) 

The difference between Ao = and A = 1 — J2& models manifests only for very small k (Fig.l) which 
gives only negligible contribution to the resulting ST/T below. Thus we use the same fits in both 
cases. 

In our present models, there are four free parameters, i.e., the primordial power-law index 
of entropy ( or baryon number) fluctuations n, the baryon density parameter (= fio in the 
present models), the dimensionless cosmological constant Ao, and finally the dimensionless Hubble 
constant h. More precisely, we adopt (|S"(A;)| 2 ) oc k n , rather than oc k n , where k = \Jk 2 + K with 
K being the curvature constant. In particular, we explore the constraints on n — fib plane for 

0. 1 < fib < 1 and —2 < n < 0. As for the remaining two parameters, we consider typical four sets, 

1. e.,(A , h) = (0, 0.5), (0, 1.0), (1 - J2 6 , 0.5), and (1 - J2 6 , 1.0). 



3. Temperature anisotropies 

3.1. Method of computing temperature fluctuations 

The linear part of ST/T in the microwave background is obtained by numerical integration of 
the linearized Boltzmann equations for photons, baryons, and massless neutrinos on the basis of the 
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gauge-invariant formalism (see, e.g., Sugiyama & Gouda 1992). The temperature anisotropy Ay is 
expanded in Fourier space by writing: 



A T (r ] ^,k) = Y J (-i) t (2l 



l)at(rj,k)Pt(fi), 



(3.1) 



where 77 is conformal time, fj, is direction cosine and Pi is the Legendre polynomial of the l-th 
order. Then the Boltzmann equation is reduced to recursive equations of a,£. We do not use a 
free-streaming approximation but rather directly solve the equations until the present epoch 770- 
Hence our computation fully includes all linear terms such as the Doppler (first-order velocity) 
effect, the Sachs- Wolfe effect and the line of sight integration of the derivative of the gravitational 
potential. Therefore our remaining task is simply to add the second-order velocity-induced terms 
to the primary ST/T so computed. In doing so, we basically follow the prescription described by 
Efstathiou (1988) except that we treat the growth rate in fio < 1 universes more properly. We 
summarize the final expressions below following his notation. 

The velocity induced temperature fluctuations can be evaluated by solving the radiation transfer 
equations perturbatively. The second-order velocity induced part of the direction averaged power- 
spectrum of the temperature fluctuations W T 2 is expressed as 



Wr 



v2 



h 

h(k) 



,v2 



2j-i 

4V x P 2 (k, 77, 



no 1 
4 

OO pi 



T yrj ,fi,k)\ dfi 
h h(k), 



3.2 



[g(vo,v')vc] Vc d ( 



2 2 M v 



Vc 



, , (1 - /i 2 )(l - 2/iy) 2 P [*Vl + V 2 - 2 M2/, Vc] P(ky, 

dyd^—rr, — 2 — o — ^2 

J-i {l + y l -2nyy 



P(k, 77c 



P(k, i] c 



f3.3) 



(3.4) 



where the suffix c represents the time at z c , D(r]) is the linear growth rate of matter density 
fluctuations, and 

r ri ] 

g(r),r)') = [axned]^ axp - a T n e a dr] . (3.5) 

L Jr\' 

The previous work by Efstathiou and Bond (1987) and Efstathiou (1988) approximated the numer- 
ical factors I\ as 6.38 irrespectively of the underlying cosmology. In the following computation, we 
numerically evaluate I\ taking account of the dependence on fio and Ao- For instance, I\ becomes 
1.68 (3.97) for h = 0.5, fio = 0.1 and Ao = (0.9) . Thus the overall amplitude of the velocity 
induced 6T/T differs by a factor of 1.95 (1.27). 

As a check of our code, we have compared our numerical results with the available data computed 
by other investigators and found good agreement; with Efstathiou and Bond (1987), Peebles (1993) 
and Stompor (1993) for transfer function, and with Efstathiou and Bond (1987) for CMB power 
spectra. 
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3.2. Constraints on fi& and n 



In order to make a direct comparison between the CMB power spectrum and specific observa- 
tions, we introduce the coefficients of the CMB anisotropies in the i^-mode defined as 



The expected temperature anisotropy for each experiment is expressed in terms of Ci and the 
specific window function Wf. 



We normalize the fluctuation amplitude by means of the COBE ST/T at 10° scale (Smoot et al. 
1992): 



First of all, in order to exhibit the resulting amplitude of the baryon density fluctuations in these 
models, we display the contour plots of the root-mean-square mass fluctuation ag at r = 8/i _1 Mpc 
for a top-hat sphere (Fig. 2). As is clearly indicated in Figure 2, <jg is very sensitive to the spectral 
index n. 

Next, with the above amplitude normalization, we predict the temperature anisotropies expected 
in specific CMB experiments. We consider the COBE quadrupole (Smoot et al. 1992), OVRO 
(Readhead et al. 1989), 9 points South Pole (SP91 : Gaier et al. 1992) and 13 points South 
Pole (SP93 : Schuster et al. 1993) experiments. As for the SP91 experiments, we only consider 
the highest frequency channel because the data in the other channels might be contaminated by 
foreground galactic emission (Gaier et al. 1992). The window function Wi of each experiment and 
Ci of the model with fl^ = 0.2, h = 0.5 and n = are shown in Figure 3. The contribution of the 
second-order velocity term is also presented. Apparently, we expect large temperature fluctuations 
on the scale of SP experiments for this model. The second-order velocity term is dominated and 
only important on OVRO scale. In Figure 4 are shown the contour plots for the COBE quadrupole, 
OVRO, SP91 and SP93 experiments on n — fl^ plane. The value of the quadrupole anisotropy 
detected by COBE is (13 ± 4) /iK. Our expected quadrupole anisotropies for desirable models, i.e., 
fib = (0.1 ~ 0.3), and n = ( — 1 ~ 0) are typically smaller than this value. If we take into account 
cosmic variance, however, these expected values are still consistent with observed one. To obtain 
constraints from OVRO and SP experiments, we use Bayesian method (Bond et al. 1991) and assume 
a uniform prior probability density. The limits from the OVRO experiment is generally weaker 
than the constraints from the SP data. It is interesting to note that the completely independent 
experiments with different instruments and the effective scales of observation consistently point to 
n ;$ — 1.5 as an allowed region. The PIB models with Ao = require that fl^ £ 0.5 additionally, 
which is in serious conflict with the prediction of the standard big-bang nucleosynthesis (SBBN), 
i.e., fifc/i 2 = (0.01 ~ 0.02). With nonvanishing Ao, especially Ao = 1 — J2& as we consider here, 
the PIB models are consistent with the SBBN prediction if n ^ — 1. The problem with such 
models, however, is that they predict very low values of <jg implying that visible galaxies should 




(3.6) 




(3.7) 



AT = (30 ± 5)/j>K, — = (1.1 ± 0.2) X 10" 5 . 



(3.8) 
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be significantly biased with respect to the underlying mass distribution. Since the PIB models are 
motivated by the assumption that non-baryonic dark matter does not exist, such conclusions seem 
quite contrived even if not excluded. 

4. Conclusions and discussion 

We have computed the ST/T in the PIB models to see whether or not the recent data including 
COBE, SP, and OVRO observations can be reconciled in a consistent set of cosmological parameters. 
As we have shown in the above, the present results constrain the range of parameters in the PIB 
models very stringently, and the models should be excluded. This conclusion, however, is subject 
to change by considering the following possibilities; the first is the uncertainty due to the cosmic 
variance. This is especially important in interpreting the quadrupole amplitude from the COBE 
data. In addition, a relatively small sky coverage of the ground-based experiments including SP91 
and SP93 might result in statistical ambiguity to some extent (Bunn et al. 1993). A clear resolution 
to the problem requires future observations with independent and large sky coverage and possibly 
at multi- wavebands. Second, we assumed in the present analysis for simplicity that the universe 
never recombined. Apparently this is not realistic. This assumption results in the most significant 
erasure of the primary ST/T , but at the same time the magnitude of the velocity-induced ST/T 
becomes largest. Therefore it is not clear whether the overall prediction for the ST/T increases or 
decreases when more realistic thermal history is considered. This should be worked out in details 
with a model of reionization after the standard recombination epoch. 
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FIGURE CAPTIONS 

Figure 1 : Transfer functions of baryon density fluctuations, (a) open models (Ao = 0); (b) flat 
models (Ao = 1 — 

Figure 2 : Contour map of <jg normalized by the COBE 10° result, (a) h=0.5, Ao = 1 — (b) 

h=0.5, A = 0; (c) h=1.0, A = 1 - O b ; (d) h=1.0, A = 0. 
Figure 3 : (a)Power spectrum of temperature anisotropies 1(21+ l)Ci/4ir for the model in which 

tt b = 0.2, h = 0.5 and n = 0. (b)Window functions W e for the COBE, SP91 and OVRO 

experiments. 

Figure 4 : Constrains on n and from OVRO, SP91, SP93 and COBE quadrupole anisotropy 
normalized by the COBE 10° result, (a) h=0.5, A = 1 - O b ; (b) h=0.5, A = 0; (c) h=1.0, 
A = 1 - O b ; (d) h=1.0, A = 0. OV(thick solid lines), Sl(thick dotted lines) and S2(thick 
short-dashed lines) are the 95% confidence-limit from OVRO, SP91 and SP93 experiments, 
respectively. Regions to the right of each curve are excluded by the corresponding experi- 
mental results. Solid lines labeled as 7, 10 and 13 are the countour curves for the COBE 
quadrupole anisotropies in units of jiK. 



